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A multi-high-frequency electron paramagnetic resonance method is used to 
probe the magnetic excitations of a dimer of single-molecule magnets. The 
measured spectra display well resolved quantum transitions involving coher- 
ent superposition states of both molecules. The behavior may be understood in 
terms of an isotropic superexchange coupling between pairs of single-molecule 
magnets, in analogy with several recently proposed quantum devices based on 
artificially fabricated quantum dots or clusters. These findings highlight the 
potential utility of supramolecular chemistry in the design of future quantum 
devices based on molecular nanomagnets. 
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Considerable effort has focused on finding building blocks with which to construct the quan- 
tum logic gates (qubits) necessary for a quantum computer (|7]|2l). Most proposals utilizing elec- 
tronic spin states take advantage of nano-fabrication methods to create artificial molecules, or 
magnetic quantum dots (EHU). A Heisenberg-type exchange coupling between dots is achieved 
by allowing the electronic wavefunctions to leak from one dot to the next. It is this coupling 
which is the essential ingredient in a quantum device because, unlike classical binary logic, 
it enables encoding of data via arbitrary superpositions of pure quantum states, e.g. |0) and 
|1) ©. These superposition states can store information far more efficiently than a classical 
binary memory. Furthermore, they permit massively parallel computations, i.e. many simul- 
taneous quantum logic operations may be implemented on a single superposition state. For a 
quantum device to become a viable technology, it should be possible to perform a reasonably 
large number of quantum operations (~ 10 4 ) on a single qubit without the superposition states 
losing phase coherence. Herein lies one of the main technical challenges, as most quantum 
systems are highly susceptible to decoherence through coupling to their environment (5). 

We demonstrate that single-molecule magnets (SMMs) may be assembled to form coupled 
quantum systems of dimers (or chains, etc.), with many of the attributes of quantum-dot-based 
schemes. Most importantly, our electron paramagnetic resonance (EPR) investigations of crys- 
tals (large, highly ordered 3D arrays) containing exchange-coupled dimers of SMMs show that 
decoherence rates are considerably less than the characteristic quantum splittings (A/h ~ GHz, 
where A is the energy splitting and h is the Planck constant) induced by the exchange couplings 
within the dimers, representing a step forward in the drive towards potential applications involv- 
ing molecular magnets. Several proposals have suggested possible quantum computing schemes 
utilizing molecular magnets The supramolecular (or "bottom-up") approach to materi- 

als design is particularly attractive, as it affords control over many key parameters required for 
a viable qubit: simple basis states may be realized through the choice of molecule; exchange 
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couplings may then be selectively designed into crystalline arrays of these molecules; finally, 
one can isolate the qubits to some degree by attaching bulky organic groups to their periphery. 

The subject of this investigation is the compound [Mn 4 03Cl 4 (0 2 CEt)3(py) 3 ] 2 -2C 6 H 14 (here- 
after [Mn 4 ] 2 ; EtCO^ is propionate, py is pyridine, and C 6 H 14 is hexane) ©, a member of a 
growing family of Mn 4 complexes which act as SMMs having a well defined ground 

state spin of S = |. This compound crystallizes in a hexagonal space group (R3) with the Mn 4 
molecules lying head-to-head on a crystallographic S$ axis. The resulting [Mn 4 ] 2 supramolec- 
ular dimer is held together by six weak C— H- • -CI hydrogen bonds (Fig. 1), leading to an 
appreciable antiferromagnetic superexchange coupling ( J ~ 10 yueV) between the Mn 4 units 
within the dimer, which influences the low-temperature quantum properties of related [Mn 4 ] 2 
dimers (1721) . Like all SMMs, [Mn 4 ] 2 displays superparamagnetic-like behavior at high tem- 
peratures, and magnetic hysteresis below a characteristic blocking temperature (~ 1 K). The 
hysteresis loops exhibit steps, which are due to magnetic quantum tunneling (MQT). However, 
unlike isolated SMMs, there is an absence of MQT at zero-field, due to a static exchange bias 
field which each molecule experiences due to its neighbor within the dimer (|72b . The effect of 
the bias is to shift the field positions of the main MQT steps by an amount of order —JS 2 / ll 
(where li is the magnetic moment of a Mn 4 monomer), so that the first step is observed on 
the hysteresis loop before reaching zero-field. However, the exchange bias by itself does not 
quantum mechanically couple the SMMs within the dimer. 

Before presenting experimental evidence for the coupled nature of the dimers, we develop 
a quantum mechanical model which takes this coupling into account. Neglecting off-diagonal 
crystal field terms and inter-molecular interactions, the effective spin Hamiltonian (to fourth 
order) for a magnetic field (B z ) applied parallel to the easy (z-) axis of a single isolated SMM 
has the form (|77J) 
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Hi = DS 2 zl + BlOl + g z fi B B z S zi , (1) 

where S zi is the z-axis spin projection operator, and the index i (— 1, 2) is used to label the two 
Mn 4 molecules in the dimer for the interacting case below; D (< 0) is the uniaxial anisotropy 
constant; -B4O4 characterizes the fourth order axial anisotropy; and g z is the z-component of the 
Lande g-tensor. The omission of transverse terms in Eq.[l]does not affect the EPR spectra (they 
merely result in weak avoided level crossings which cause the MQT). 

For the case of two quantum mechanically coupled SMMs, the effective dimer Hamiltonian 
(Hd) may be separated into the following diagonal and off-diagonal terms 

H D = [4 + H 2 + J z S zl S z2 ] + {\J xy {StS 2 + S^S+)}, (2) 

where H\ and H 2 are given by Eq. [H the cross terms describe the exchange coupling between 
the two SMMs within the dimer, and the J values characterize the strength of this coupling. 
The diagonal zeroth order Hamiltonian (H 0D , in square brackets) includes the exchange bias 
J z S z \S z2 which has been considered previously (1721 . The zeroth order eigenvectors for the 
dimer may be written as products of the single-molecule eigenvectors \mi) and \m 2 ) (abbrevi- 
ated \mi,m 2 )), where m\ and m 2 represent the spin projections of the two molecules within 
the dimer. The zeroth order eigenvalues are then easily obtained by solving Eq. 1 separately for 
molecules 1 and 2, and adding the exchange bias J z m 1 m 2 (Fig. 2). 

In EPR, the only effect of the exchange bias is to cause shifts in the positions (energies) 
of single-spin transitions (mj — > rrii ± 1), with the magnitude of the shift (bias) depending 
on the state mj of the other molecule within the dimer. It is the off-diagonal interaction in 
Eq. 13 (H' D in curly brackets) which couples the molecules, giving rise to the possibility of 
single-photon transitions between coupled states of the dimer. In principle, one could observe 
this coupling in hysteresis measurements, as magnetic relaxation mediated by tunneling into 
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distinct superposition states occurs at slightly different magnetic field strengths, even when the 
tunneling occurs via states with the same total spin projection (= mi + m 2 ). However, the 
predicted splittings of MQT resonances turn out to be less than the inhomogeneous linewidths 
of the hysteresis steps (1721) . Thus, clear evidence for the coupled nature of the dimer system has 
so far been lacking. 

For illustrative purposes, we treat H' D perturbatively. As this interaction conserves angular 
momentum, the eigenvectors may be grouped into multiplets based on the sum of the projec- 
tions M = mi + m 2 . H' D then acts only between states within a given multiplet. The zeroth 
order eigenvectors are grouped according to this scheme in Fig 2 (left). In first order, H' D acts 
between zeroth order eigenvectors |mi, m 2 ) and \m\ ± 1, m 2 =F 1). The effect of this first order 
interaction is most apparent in the M = — 8 multiplet, where it lifts the degeneracy between 
the | — |, — |) and | — |, — |) states. The resultant eigenvectors correspond to symmetric (S) 
and antisymmetric (A) superpositions of the original product states. Indeed, H' D causes consid- 
erable mixing of the zeroth order eigenvectors within all multiplets, resulting in the first order 
corrected eigenvectors which are listed in Fig. 2 (right) for the lowest four multiplets; here, 
l m i! m 2)s implies (|mi, m 2 ) + |m 2 , mi)) and \m\, m 2 )A implies (|mi, m 2 ) — |m 2 , mi)). 

In Fig. 2, we display a schematic of the energy level shifts and splittings (not to scale) caused 
by the exchange bias, and by the full exchange, for the lowest lying levels at high magnetic fields 
(M = —9 to —6). The states are numbered for convenient discussion of the data. For clarity, 
higher lying states with M > —6, including the zero-field | ± |, =f|) ground states, are not 
shown in Fig. 2. Application of a magnetic field parallel to the easy axis merely shifts all of the 
zeroth order levels by an amount gfi B B z M. Thus, 5M= ±1 EPR transition matrix elements may 
be accurately calculated using the eigenvectors in Fig. 2. The magnetic dipole perturbation only 
allows transitions between states having the same symmetry. The strongest of these transitions 
are shown in Fig. 2, labeled (a) through (g). 
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In the left-hand panel of Fig. 3, we display temperature dependent high-frequency EPR 
spectra obtained at 145 GHz, with the magnetic field applied parallel to the easy (z-) axis of 
a small (< 1 mm 3 ) single-crystal sample; details concerning our high-frequency EPR setup 
are given elsewhere (T7H). The inset shows a single 6 K spectrum (/ = 140 GHz) for a related 
monomeric Mn 4 complex without head-to-head interactions (|7?|) . The monomer data are typical 
of most SMMs, showing a series of more-or-less evenly spaced resonances, and a smooth vari- 
ation in intensity from one peak to the next. By contrast, the dimer spectra exhibit considerable 
complexity. In spite of this, the simulated dimer spectra (colored traces in the right-hand panel 
of Fig. 3) show remarkable agreement with the raw data, both in the peak positions and relative 
intensities. The optimum parameters were deduced from a single fit to Eq.|2]of the main EPR 
peak positions obtained at many microwave frequencies. This fit, displayed in Fig. 4, yields the 
following values: D = -0.750(15) K, 5° = -5(2) x 10~ 5 K, g z = 2 and J = 0.12(1) K. 
These crystal field parameters are very similar to those obtained for the monomer [D = —0.7 K, 
B® = — 9 x 10~ 5 K (I7?l) l. We did not find it necessary to include anisotropy in the superex- 
change coupling for the dimer (i.e. J z = J xy = J), though long range dipolar interactions 
improved the quality of the fit <\15\\16h . 

The simulated spectra (Fig. 3) are mainly limited to transitions among the levels displayed 
in Fig. 2 [(a) through (g)]; we have also included the (7)s,a — > | — |, — \) and | — |, — |) — > 
| — |, +|) transitions, labeled (h) and (i) respectively. Resonance (x), meanwhile, corresponds 
to the degenerate | + §,— §) — > | + §,— |) and | — §,+§) — > | — |>+|) transitions. The 
only significant differences between the experimental data and simulated spectra are seen in the 
2 — 3 T region, which is due to fact that we did not consider several moderately strong transitions 
involving higher lying (M > —6) states. We deliberately avoid reference to superposition states 
in discussing resonance (x), as the interaction between the | ± |, =f|) states is extremely weak 
(9th order in H' D ). Consequently, even the weakest coupling to the environment would likely 
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destroy any coherence associated with the 2 _1 / 2 | + |, —%)s,a superposition states. 

Resonance (x) is observed only over a narrow low-field region (< 0.7 T) over which the 
I ^ §> "F§ ) levels represent the ground states of the dimer. By following the relative intensities 
of resonances (x) and (a), one obtains an independent thermodynamic estimate of the exchange 
bias which is in excellent agreement with the value obtained above, and with independent hys- 
teresis measurements for the same complex (17). We note that the previously published mea- 
surements of the exchange bias in [Mn 4 ] 2 involved a slightly different solvent of crystallization, 
the full compound having the form [M^OgCl^OaCEtMpy^-SMeCN (|P1I23- EPR studies 
for this complex (14) show fewer transitions from excited levels [transitions (x), (a), (b) and (c) 
remain clearly visible]. Nevertheless, one can still estimate a coupling constant [J = 0.10(1) K] 
from the exchange bias, which is in agreement with the published value (U2l . 

The inset to the right panel of Fig. 3 shows that it is the transverse part of the exchange {H' D ) 
which brings the simulations into excellent agreement with the data. Indeed, there is no way to 
obtain anything closely resembling the experimental data without including H' D in the calcula- 
tion, thus providing compelling evidence that the molecules are coupled quantum mechanically. 
The issue of quantum coherence is best illustrated by examining the splitting of resonances (f) 
and (g) — this splitting is directly proportional to J xy , and corresponds to the ~ 9 GHz shift 
of the (4)5 level relative to (5)^ (Fig. 2). If the phase decoherence rate (r^ 1 = characteristic 
rate associated with the collapse of a quantum mechanical superposition state) were to exceed 
9 GHz, one would expect broad EPR peaks due to transitions between bands of incoherent 
states; these bands would occupy the gaps between the energies given by the exchange bias pic- 
ture and the full exchange calculation in Fig. 2, thereby smearing out most of the sharp features 
in the observed spectrum. In principle, is the same as the transverse spin relaxation time T 2 , 
which can be estimated from EPR linewidths (AM = ±1 transitions) (18). However, we know 
that these widths are dominated by weak dimer-to-dimer variations (strains) in the Hamiltonian 
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parameters, i.e. the actual r^ 1 is buried within the inhomogeneous EPR linewidths d7?lli5lli(5l) . 
and is probably much less than 9 GHz. As a worst case, the narrowest EPR lines would imply a 
decoherence time on the order of 1 ns. In order to determine the real T 2 (= r^), one should carry 
out time resolved (pulsed) EPR experiments, e.g. the free-induction-decay of an initially satu- 
rated EPR transition, or Rabi spectroscopy dTffl) . Time resolved experiments in this frequency 
range are technically challenging but, nevertheless, represent a future objective. 

The magnitudes of the quantum splittings (in frequency units) provide a rough estimate of 
the rates at which one could perform computations. In comparison to many competing technolo- 
gies [e.g. NMR (19)] these rates are high for electronic spin states, i.e. GHz rather than kHz or 
MHz. The largest quantum splittings (A/h) for the dimer are on the order of a few tens of GHz. 
In fact, Ar^/h represents a rough figure of merit for a quantum device, as it gives an estimate 
of the number of qubit operations one could perform without loss of phase coherence. For the 
worst case given above, Ar^/h ~ 30 — 100; in reality, it may well be 10 4 or greater. The most 
useful coupled states of the dimer would be the antiferromagnet zero-field 2~ l l 2 \ + |, — §)s,a 
ground states, or Bell states (2). As already discussed, the tunnel splitting of these states is 
negligible in zero-field (~ Hz). However, it is possible to increase this splitting to a practical 
range (~ GHz) with a transverse magnetic field. While there remain technical challenges along 
the road map towards molecule-based quantum devices (e.g. low operating temperatures, meth- 
ods for addressing nanometer-sized molecules, etc.), the present study demonstrates that the 
"bottom-up" (molecular) approach provides excellent opportunities to study coherent quantum 
superposition states. Future materials design strategies will, therefore, explore the following 
possibilities: optical control of the exchange coupling between the two halves of a dimer; in- 
creased isolation of the dimers in order to further reduce decoherence; and the inclusion of some 
form of asymmetry within the dimer (e.g. uncompensated electronic spins, or selective nuclear 
spin labeling), thereby facilitating readout of the state of the system. 
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Fig. 1. The [MiuOsCUtC^CEtMpy^h dimer; the Mn 3+ and Mn 4+ ions are labeled Mn and 
Mn', respectively. The dashed lines represent the C— H- • -CI hydrogen bonds holding the dimer 
together, and the dotted line is the close approach of the central bridging CI atoms believed to 
be the main pathway for the exchange interaction between the two Mn 4 molecules. 

Fig. 2. Schematic showing the lowest energy states (M= —9 to —6, not to scale) of the dimer: 
the zeroth order energy levels and eigenvectors are shown on the left; energy shifts due to the 
exchange bias are shown in the center; and the results of a full quantum calculation (Eq. |3) are 
displayed on the right. The colors denote the total angular momentum (M) state of the dimer, 
and the levels have been numbered to aid discussion. Several of the strongest EPR transitions 
are indicated by arrows [(a) thru (g)]. The corrected eigenvectors are listed next to each state: 
based on the deduced value of J {vide infra), (3 = 0.400 and j3' = 0.231; the subscripts S 
and A (on a state |rai, m 2 )) respectively denote symmetric and antisymmetric combinations of 
\mi, m 2 ) and |m 2 , mi); and the C s are normalization constants. 

Fig. 3. The left-hand panel displays temperature dependent easy-axis data obtained for the 
[Mn 4 ] 2 dimer at 145 GHz (the dips in transmission correspond to EPR); the inset (black trace) 
shows a single 6 K, 140 GHz, spectrum obtained for a monomeric Mn 4 complex [the resonances 
are labeled according to the m states from which the transitions were excited (|7?)) 1. The right- 
hand panel contains simulations of the dimer data, while the inset illustrates the effect of the 
transverse part of the exchange (J xy ) for four values of J xy /J z (T= 8 K). In both figures, 
resonances (a) through (g) correspond to the labeled transitions in Fig. 2; resonances (x), (h) 
and (i) are discussed in the main text. A Gaussian distribution in D (a D ~ 1%) was included in 
the simulations in order to obtain realistic lineshapes (1751). 

Fig. 4. A single fit to Eq.[2]of the positions of EPR peaks obtained at several frequencies. The 
optimum Hamiltonian parameters were obtained from this fit. The transitions have been labeled 
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for comparisons with Figs. 2 and 3. 
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Figure 2: S. Hill et al. 
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